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Development of the Disturbing Function in Planetary Theory , in 
terms of the mean anomalies and constant elliptic elements. 

By P. H. Cowell, M.A., F.R.S. 

1. The following notation will be employed :— 

a mean distance 
/ true anomaly 
r radius vector 
g mean anomaly 

e = tan \ ifr where sin i j/ is the eccentricity. 

Plain letters refer to the disturbed planet, accented letters to 
the disturbing planet. 

o), a) the distance from the intersection of .the orbits to the 
perihelia. 

o- = sine of half the inclination of the two orbit planes. 

2, Then the disturbing function (divided by the mass of the 
disturbing planet) is 

R = {r 2 + r' 2 - 2 rr cos S} ~& - 7 CQ ^ ^ 

where cos S = cos ( -/+/' - w + w) - <r 2 {cos ( -/+ f' ~ <o + (o') - 

cos (/+/+ oj + w')} 

If in this expression we replace r by a, f by a\ we have 
R = {a? + a 2 - 2aa' cos S} 

n. & 


3. Suppose we have an expression 

</>(a)ePV- 1 

and we require to replace a by r and / by its value in terms of g , 
and we wish to have the modified function expressed in a series of 
multiples of we may adopt the following notation : 


-rw d 

D = a— 
da 


j 


1 = 


22 




{n t j, I, D}g( I +.?)0v'- 1 


where the quantity in brackets {n,j, I, D} is an operator operating 
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on <£(a) and the summation extends over all values of n> j 
that make the quantities under the factorial sign zero or positive 
integers. 

Now suppose we have an expression 

<£(a, 

for the expansion of 

4 >(r, r')e^ + 1 'f^ ~ L ; 

we then obtain 






'ri 




{n.jl.V ).} { n'.j .I'.D'.}«[(! +j) 9 +a'+?)g'W 


the summation again extending over all values of n, j , n, j' that 
make the quantities under the factorial signs zero or positive 
integers. 

The advantage of employing the purely numerical denominator 
shown above is that the coefficients in the expanded values of 
{w, j , I, x} are integers. 

These operators are, moreover, rational integral algebraic 
functions of I, x of order n. 

4. The values of {n, j, I, x }, as far as n = 6 , are given in the 
following table. 

We have {n. - j . I. x.} — {n.j. - I. x}. 


{o, o, I, x} = 1 

{1, 1, I, x} = 4 - 2I - x 

{2, 2, I, x) = +5l + 4l 2 -3a;-4l^ + ^ 2 

{2, o, I, sc} = - 4I 2 4 - x +x 2 

5. These formulae may be very easily developed from special 
cases. For instance, put x=o, 1 = 1, then we have 

€ n 

00s I, 0} cos (g +jg + 6 ) 


or 




n+j\ n-j 


{w, j\ 1, oJeb+.hsV- 
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where 
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If this be regarded as a known elliptic formula, to be used as a 
starting-point, by raising to the I th power we have 


^"“22 



{w, j, I, o}eP~ + ^9V-T- 


thus verifying all the special cases of the table obtained by putting 
x = o. 

6. Again, put I = o and we have 

= 2 2 ^+^* °' D} ^ (a) G08jg 

where 


{0 

O 

0 


I 






{1 

I 

0 
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D} 
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3D + 

D 2 
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D + 

D 2 




{3 

3 

0 

D} = 
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0 

D} 

{4 

0 

0 

D} = 

- 

6D- 

9D 2 - 

2D 3 + D 4 



Is 

5 

0 

D} = 

- 

1569D + 

1220D 2 - 

305D 3 + 30D 4 - D 6 -{5 -5 

0 

D} 

{5 

3 

0 

D} = 

+ 

435 D- 

176D 2 - 

29D 3 + 14D 4 - D 6 = {5 — 3 

0 

D} 

{5 

1 

0 

D} = 

- 

63D- 

26D 2 + 

i3D 8 + 6D 4 - D 6 ={5-1 

0 

D} 

16 

6 

0 

D} = 

- 

21576D +18694D 2 - 

5595 l )3 + 745 l )4 - 45 l )6 + D 6 ={6 -6 

0 

D} 

{6 

4 

0 

D} = 

f 5936D- 

3366D 2 + 

io5D 3 +165D 4 -25D 6 + D 6 = {6 -4 

0 

D} 

{6 

2 

0 

D} = 

- 

992D + 

46D 2 + 

I97D 3 + 9D 4 -i3D 5 + D 6 ={6 -2 

0 

D } 

{6 

0 

0 

D} = 

+ 

60D + 

1 l8D 2 + 

57D 3 - 11D 4 - 9D 6 + D 6 




© Royal Astronomical Society • Provided by the NASA Astrophysics Data System 


Downloaded from http://mnras.oxfordjournals.org/ by guest on March 24, 2015 









1909MNRAS..69..170C 


174 Mr. P. H. Cowell, Development of the lxix. 3, 

but if log r=log a + p, 

we know, by Maclaurin’s Theorem, that 

<t>(r) = 2 ^ D to 
m ! 

Hence, in the above set of formulae, picking out the coefficient of 
the first power of D, we have, on combining cos jg with cos -jg, 


p= “ 

e 

cos g 

x 2 

- 

-1 e 2 

3 cos 2 g 

x 2 

4 - 

€ 2 



- 


17 cos $g 

x 2 

+ 

*«* 

5 cosp 

x 2 

- 

1 *4 
2T 6 

142 cos 4 g 

x 2 

+ 

1 f 4 

40 cos 2/7 

x 2 

- 




- 

1 f-b 

T2 (7 6 

1569 cos 5^ 

x 2 

+ 

1 

2T 6 

435 cos 3 ^ 

x 2 

- 

_A_ e 5 

12 6 

63 cos g 

x 2 


- e 6 21576 cos 6# x 2 
-bnhr* 6 5936 cos 4^x2 

- e 6 992 cos 2g x 2 


If this is regarded as a known elliptic formula we may calculate 


for m= 2, 3, 4, etc. 


and thus verify the values of 


{n. j. o. D} 


7 * • 

7. Taking values of cos / and log — in terms of c and g as 

a 

known, we have now shown how to obtain 

{n,j, o, D} and {n,j, I, o} 

We have lastly to obtain 

{»,/> I, D} 

We are supposed to start from <£(&) eW- 1 in which we 
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require to replace a by r and / by its value in terms of g, and we 
then wish to pick out the coefficient of 


n +A \ ( U ~J 


_g(I-h l)9V ~ j 


We have seen how to replace </>(a) by a series whose typical 
term is 


{v. y.o.DjeYffS- 1 


v + 7 \ ( v ~7 


and eb/V-i by a series whose typical term is 


n-v+ 3 -y \ ( n-v-j+y \ 


{n — v,j— y, I. o}e( I+ i~y') 9^ / ~ 1 


Hence by multiplication we see that 

'n+j\\ ( n ~j 

,^ 1 = 22/7 


±y\\ fizy\( n ~ v +J~y\ (n-r-J+y 


'{v.y.o.D }{n- 


the summation extending over all values of v, y which make 
every quantity under the factorial sign in the denominator zero 
or a positive integer. 

8 . We have now merely to throw 

E ,= {a 2 + a 2 - 2 aa! cos S} ■£- —^ ^ 


into the form 

2 </>(«, «') cos {lf+Vf + 6 } 

where 0 , a combination of <d and a/, is a constant. 

The second term of R is already in the required form, directly 
the value of cos S is substituted from § 2 of this paper. 

We also have 

{a 2 + a ' 2 - 2 ad cos S} 

= 2 { « 2 + «' 2 ~ 2aa ' cos (-/ +/ - «> + <■>')} (aa'f 

Jc =O 2 W 


x cr 2k {cos (/+/' +co + <o) - cos ( -f+f - ft)+ (!)')}* 


v-j- 
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Now if a be the larger of the two quantities a, a', and if a be the 
ratio of the smaller to the larger, 

a {a? + d 2 - 2 ad cos 0}- k ~^{aa!) k +a k 

= { 1 — ae 0, v / -i| - — ae~ 9 ^~ 1 ]■-***-4 

- (, +'*±1 J- / -‘ + ( *+*>■(*+ . 


I . 2 


- 9 \/ - 1 


X 


ae 


+ 


(* +*)•(*+ kW * 


(■+¥ 


= i+(t+i)\* + (*+m*+i) a n * + .... 


+ 


I . 2 


I 2 . 2 2 


+ ^ + ^ + (^4-1).(&4-f) ff2 + {(^ + i)(^ + #)} , {(fc + l)(^ + f)} a 4 + 

2 -3 


I . 2 

+ etc. 


1 . 2 


cos 0 + cos 


= 2, &i.fc COS a . 0 where b-i.k = h.k 

% — CO 

We now easily get 

1 — 03 K — 00 

r= 2,2°' 2 ' cA «-' . C0s{t( —f +f — u) + g/) + /c(/ + c /’ / + 0 ) + a/)} 


1 = - 00 /c = o 


where, sufficiently for a development as far as or 6 , 
cl _ 0 = a A -i. o~b%. o~ ^cP • a { b{~ x. x + &£-fi. 1} 

+ f cr 4 a 2 {^-2.2 + 4 ^’. 2 + bi + 2 . 2 } 

~ -^Or 6 a 3 {^_3,3 + I • 3 4 " 9 &t+i. 3 4 * ^i+3 . 3} 
aAi . X = CL b% # ! “ 7 T (Po?(bi _ X . 2 4 " I « 2) 

4 * _ 2 '. 3 4 " 3^. 3 4 " ^i-j-2.3) 

&Ai. 2 = \ct?bi _ 2 -g"< Po?{bi _ x , 3 4 ~ ^i+i. 3) 

aAi. 3 = fa 3 ^. 3 . 

On account of the second term of R we must add 


1 ci 

--- 72(1 -or 2 ) to Ax.o and A_ x . 0 

2 iv 


and 


a 


a 


'2 


to Ao. x. 
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9. Summary :— 

The value of the disturbing function can therefore be expressed 
in the perfectly general form 

%—ca k= co 

^ ^ 2 n ' 2/ 

i= - 00 Jc = o 



^n ^n' 




{n,j, -i + k, D}{«',/, + i + k, D'JA,.,. 


x cos {( — i + k 4 - (i + k-\-j + ( — i 4 - k)w + ( i + &)<*/} 

where since R is of the order - i with respect to a, a' we may 
replace D' by - i - D if we please, and where the value of every 
quantity required for an expansion as far as € 6 , e' 6 , cr 6 (including 
such combinations as eW of the eighteenth order) is given in this 
paper. 

10. Numerical example :— * 

Required the coefficient of 

2f r! i T r i cos (" 29 + 49 ~ 3 “ + 5"') 

Here 

i = 4 k= 1 w = 3 j = 1 n = 1 j'— - 1 
The two operators are 

{3.1. - 3 . D} = + 144 + 26D - 5D 2 - D 8 
{1, -I, +S» D '} = { 1 . + 1 . - 5 > D '}= - io-D'= -9 + D 
The quantity operated upon is o- 2 A 4 . l 

= — | a & 4 ,1 - — o- 2 a 2 (6 3 2 + b 51 ) l 

a 1 2 j 


stopping at cr 4 c 3 c 1 or the eighth order 



6. 7. 9 

• 2 * 2 * ~ 5 T 

. 2 . 3.4 



1JL 3 11 13 3 5 

+ + 2 J_ f a 4 + 

5 1 5.61.2 


_6 7 _9 

A — 2 * ^ • 2„3 

O3.2- a 


1.2.3 


5 . i 
h — 2 * 2 
^ 6.2 - 


5 . 

2 • 


1.2.3 


11 
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I+- 2 -. 

- 2 -a*+ . 
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3-3 ( 

1 6 



I + X. 

_2_a 2 + 


•4-5 1 61 
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Further, if Venus be the disturbed planet and the Earth the 
disturbing planet, then 

1 

a = 0‘723 331 5 = a 
a=i. = a 

logo-= 8*471 368 

(0 = 234° f 22 ,,# 2 • €='003 421 70 

o> / = 2o5° 1' 29 ,, *i e^‘008 386 14 

(see Newcomb, A.P. , iii. p. 406, 459) 

and 

Da w = ma m . 

The principal term becomes 


o > 


3 

¥ 


5 


X JL 
2 • 2 


— — . 96 . o-V\ 

2 I • 2 . 3.4 


which is less than 10 n , and has still to be multiplied by the mass 
of the Earth. 


Note on Major MacMahon's paper “ On the Determination of the 
apparent Diameter of a fixed Star ” By A. S. Eddington, 
M.A., M.Sc. 

In Major MacMahon's paper (p. 126 ante), the fading of the 
light of a star at the instant of disappearance behind the Moon was 
treated purely as a problem of geometrical optics. When account 
is taken of the diffraction of the star's light by the limb of the 
Moon, his results will be found to require some modification. The 
argument given in this note appears to show that, judged by the 
rate of diminution of its light at occultation, every star will seem 
to have an angular diameter of about o"*oo8, due to this 
diffraction. 

Reducing the problem to its ideally simple form, we consider 
a point-star S. The Moon is treated as a screen M M' with a 
straight edge through M perpendicular to the plane of the 
paper. 

Let X E X' be a screen in the position of the observer, E being 
the boundary of the geometrical shadow cast by the Moon. 
Let J be the intensity of illumination of the screen by the star 
when the Moon is out of the way, the true illumination of the 
screen will then be as follows: * — 

Inside the geometrical shadow, the illumination is zero until 
near the edge; then it slowly increases, reaching the value \ J at E; 
outside the geometrical shadow it still increases, until it reaches 

* See, for instance, A. Schuster, Theory of Optics, p. 94. 
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